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Abstract. A multidimensional Brownian motion with partial reflection on a hyperplane 
S in the direction qN + a, where A'^ is the conormal vector to the hyperplane and q S 
[—1,1], a g S are given parametres, is constructed and this construction is based on both 
analytic and probabilistic approaches. The joint distribution of d— dimensional analogy to 
skew Brownian motion and its local time on the hyperplane is obtained. 



§1. Introduction. 

In this paper we construct a generalized diffusion process in a d-dimensional Euclidean 
space 3?'' for which the diffusion matrix B is constant and the drift vector is equal to 
{qBi/ + a)6s{x) where i/ G 3?'' is a given unit vector, S is the hyperplane in 3?'' orthogonal 
to v, the parameters q G [^Ij 1] a-nd a €z S are given and 6s{x) is a generalized function 
that is determined by the relation 



valid for any test function ip. 

This problem in a more general case (namely, the parameters a, q are some functions 
of X e S* and the diffusion matrix i? is a function of x € K'') was considered by Kopytko 
B.I. [1] and solved by an analytic method. More simple case {B is an identical operator, 
a,q are some constants) was considered by Kopytko B.I., Portenko N.I. [2] with use of 
an analytic method too. 

In this work a probabilistic approach to this problem will be proposed. We construct 
the process desired as a solution of the following stochastic differential equation 

dx{t) = {qBv + a)6s{x{t))dt + B^/^dw{t) 



where _B^/^ is the positive square root from the operator B and w{t) is a given Wiener 
process in 3?''. If a = then the solution of this equation is known. Let us denote it by 
x{t). The idea is to add the process a Ss{x{T))dT to the projection of x{t) on S. 

The paper consists of two sections. In Section 2 we describe the analytic method. It 
is similar to that of [2]. The probabilistic approach is presented in Section 3. 
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§2. The Analytic Method. 

Let v he a fixed vector with \\u\\ = 1 and S be the hyperplane in Ji'' orthogonal 
to v 

S = {x €^'^\{x,iy) =0} 

And let x{t) be a continuous Markov process in with transition probability density 
9o{t,x,y) 

go{t,x,y) = [det{27rtB)]-'/\xp^-^^{B-\y - x),y - x)Yt > 0,x e^'',y e^'' 

where B is the given positive definite symmetric operator in 

Let g G [— 1, 1], a G S' be given. We are looking for a function u{t,x,(fi) defined for 

t> 0,x e^'^,ip e Cbi^'^) such that 

(1) u(t, X, (fi) satisfies the heat equation in the domain t > 0,x ^ S 

If we fix a coordinates in JR'' then D'^u = {g§-^)ij=i and this equation can be 
written as follows 

du 1 , d'^u 



at 



2 12 ^'"^ dxidx-i 



where B= 



(2) u{t, x, ip) satisfies the initial condition 

limu{t,x,ip) = ip{x),x €?lt'^,tp€ Cbi^H'^) 

(3) For x G S and t > the relation 

u{t, X+, if) = u{t, X—, ip) 

holds, where u{t, x±, if) = lime4,o u{t, x ± ev, ip) 

(4) The relation 

l + qdu{t,x+,ip) 1- qdu{t,x-,ip) , ^ . „ 
-2 ON 2 aA^ + ^"'^")=°' 

is valid for a; S 5 and t > , where is gradient of u{t, x, ip), N is the conormal 
vector to 5*, = Bv. 

Let us denote by Xi,{t) the projection of x{t) on i', by tts the operator of orthogonal 
projection on S such that x^{t) — nsx{t) is the projection of x{t) on S. If v is not an 
eigen vector of the operator B then x^{t) and x^{t) are not independent processes and 
we can not represent go{t,x,y) as a product of Xv{t) and x^{t) densities. But we may 
consider the following decomposition of go{t, x, y) 

/, X 1 / {Vv-x^f 
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X [det{2-KtBs)]-"^ expi^-^^{B-,\ys - xs - ^^6), ys - xs - ^^6)} 

where b = ttsBv, Bs = {ttsB~^tts)~^ ,(J^ = {By, v). 

The first factor of this decomposition is the density of the process x^{t), the second one 
is the conditional density of the process x^{t) under the condition Xi,{t) = yv Denoting 
by g^{t, X, y) the second factor we may write 

5o(t, X, y) = -^L= exp /(i, x, y), t > 0, x e JR^ y e SR"^ 

The problem (1) — (4) may be solved in the following manner. Let us fix a system of 
coordinates so that v = and make use of the Fourier transformation with respect to 
a;i,...,Xd_i and the Laplace transformation with respect to t to the function u{t,x,(p). 
It will reduce this problem to solving a pair of second order differential equations with 
constant coefficients. Then free coefficients may be found from boundary conditions 
(3), (4). Returning to the original function we get the following assertion. 

Theorem 1. If \q\ < 1 then there exists a continuous Markov process in with tran- 
sition probability density G{t, x, y) such that 

u{t, x,^)= [ ^{y)G{t, X, y)dy, i > 0, x G 3?^ ^ G Cb{3i'') 



where 



G{t,x,y) 



1 



exp i — ^ - exp 



2ta2 j ^^^^i 2ta^ 
^2/1 I u I _L i«, I r i'^2/) I u I _L i„, l^2■ 



9''{t,x,y)+ 



Remark. The condition |g| < 1 provides that G{t,x,y) is positive. 

§3. The Probabilistic Method. 

Let w{t) be a given rf-dimensional Wiener process, -B be a symmetric positive definite 
operator in 3?'* and parameters q G [—1,1], a G 5 be given. Consider the stochastic 
differential equation 

(1) dx{t) = [qBv + a)5s{x{t))dt + B^/'^dwit) 

where ds{x) is the generalized function on K'' defined above. 

When a = then the solution of the equation (1) is known. It may be called a 
multidimensional analogy to skew Brownian motion. Let us denote it by x{t). The 
density of x{t) equals to 



g{t,x,y) = 



exp|-^^|+,signy.exp| 



9 {t,x,y) 



where t > 0,x G ^'^,y G ^'^ 
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The equation (1) means that we have to find the process x{t) that satisfies the foUowing 
properties. The projection of x{t) on 1/ is skew Brownian motion, the projection of x{t) 
on S is the Gaussian process in S (with mean + ^"'2" b and correlation operator tBs) 
plus the process arjt, where rjt is the functional of the process x{t) such that 



Ss{x{T))dT 



Remark. Actually rit depends only on Xi,{t). Therefore we may write rjt = /q Ss{x{T))dT 
The functional r]t is a nonnegative continuous homogeneous additive functional, it 

increases at those instants of time for which the process x{t) hits the hyperplane S. 
To solve the equation (1) we have to find the joint distribution of x{t) and r]t. 

Lemma 1. The joint distribution of x{t) andrjt has the form 



Fx {i{t) e dy,f]t ede} = { 



s{0) 



exp 



{yu - x^Y 
2ta'^ 



exp 



{\y.\ + K\f 



+{l + q sign y^) ^^^^3^2 \ " 2ta^ — j ^ ' ^' ^^'^^'^^ 

Remark. {x{t) G dy, rit G d6} has the atom at the point ^ = 0. 
Proof. Consider 

u{t, X, A, /x) = E^e'(*(*)''')+*^''* ,t>0,xe^^,Xe^,fj,G^^ 

where E^; is the averaging operator corresponding g{t,x,y). 

The function u{t,x,X,n) satisfies the following integral equation (see [3]) 

(2) u{t,x,\,iJ,)= / e^'^^'^^g{t,x,y)dy + iX dr u{t — t, z,X, fj,)g{T,x, z)dz 

Jo Js 

Note that the function u{t, ■, X, jj) on the right hand side of (2) depends on its values 
on S. Then using Fourier-Laplace transformation we can find u{t, x, X, /z) for a; e S 

u{t,x,X,fx)= [ de [ e^(^''^)+^^'(l + gsign2/,)x 
Jo J?st<' 



V27rto2 



exp< - 



/ ia^e+\x,\ + \y,\)^ \ su ^ , 



Substituting this relation in (2) we will obtain 



u{t,x, X,iJ.) 



sjt<i-\/27rto^ 

de 



exp 



'iVu - x^f 
2ta^ 



exp 



g^{t,x,y)dy+ 



X exp 



{a^e+\x,\ + \y,\)' 

2iCT2 



a^9+\x^\ + \y^\ 
V27ri3cr2 

g^{t,x,y)dy 



where t > 0, a; e 3?"^, A e SR, e 3?'*. This gives us {x{t) e dy, r]t e dO}. 

If we know the joint distribution of x{t) and r]t, then the equation (1) can be solved. 
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Theorem 2. The soluUon of the equation (1) is a continuous Markov process x{t) with 
transition probability density G{t,x,y). 



Proof. We have 

(3) P. 



{x{t) edy}= / Px {x{t) e dy, rjt £ dO} 

^[0,+oo) 



When 6 = then x{t) = x{t) and joint distribution P^; {x{t) G dy, rjt € dO} is equal to 
Px {x{t) e dy, ijt e d6} at the point ^ = 0. For each 9 e (0, oo) the process x{t) is equal 
to x{t) + aO. Integrating (3) on 6 we obtain that 



Px {x{t) e dy} 



1 



■ / {yu - x^f \ / 



+ 



/ (^signy;, + 1)- 



■ exp 



It coincides with the result of the Theorem 1. 



2^2 



g^{t,x,y)dy 



g^{t,x + a9,y)d0dy 
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